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Abstract 

Given n points in a circular region C in the plane, we study the problems of moving the 
n points to its boundary to form a regular n-gon such that the maximum (min-max) or the 
sum (min-sum) of the Euclidean distances traveled by the points is minimized. The problems 
have applications, e.g., in mobile sensor barrier coverage of wireless sensor networks. The min- 
max problem further has two versions: the decision version and optimization version. For the 
min-max problem, we present an 0{n\og n) time algorithm for the decision version and an 

f^ I 0(nlog n) time algorithm for the optimization version. The previously best algorithms for the 

P\J . two problem versions take 0{n^'^) time and 0{n^'^ logn) time, respectively. For the min-sum 

problem, we show that a special case with all points initially lying on the boundary of the 
circular region can be solved in 0{n^) time, improving a previous 0{n^) time solution. For the 
general min-sum problem, we present a 3-approximation 0{n^) time algorithm, improving the 

^^ ■ previous (1 -I- 7r)-approximation 0{n^) time algorithm. A by-product of our techniques is an 

algorithm for dynamically maintaining the maximum matching of a circular convex bipartite 

r^ . graph; our algorithm can handle each vertex insertion or deletion on the graph in 0(log^ n) 

r \ ' time. This result is interesting in its own right. 

1 Introduction 

►^ . Given n points in a circular region C in the plane, we study the problems of moving the n points 

CNj I to its boundary to form a regular n-gon such that the maximum (min-max) or the sum (min-sum) 

of the Euclidean distances traveled by the points is minimized. The problems have applications, 
e.g., in mobile sensor barrier coverage of wireless sensor networks. The problems have been studied 
r^ , before. In this paper we present new algorithms that significantly improve the previous solutions 

^^ I for the problems. 

1.1 Problem Definitions 

^ . Let \ab\ denote the Euclidean length of the line segment with two endpoints a and b in the plane. 

H I Let C be a circular region in the plane. Given a set of n points S = {Aq,Ai, . . . , An-i} in C (i.e., 

in its interior or on its boundary), we wish to move all sensors to n points A'q, A'^, . . . , A'^_^ on the 
boundary of C that form a regular n-gon. The min-max problem aims to minimize the maximum 
Euclidean distance traveled by all points, i.e., maxo<i<ra_i{|^i^^|}. The m.in-sum problem aims to 
minimize the sum of the Euclidean distances traveled by all points, i.e., Y17=o \^i^i\- 

Further, given a value A > 0, the decision version of the min-max problem is to determine 
whether it is possible to move all points in S to the boundary of C to form a regular n-gon such 
that the distance traveled by each point is no more than A. Indeed, let Ac be the maximum distance 
traveled by the points in an optimal solution for the min-max problem. Then, the answer to the 
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feasibility problem is "yes" if and only if Ac < A. For discrimination, we refer to the original 
min-max problem as the optimization version of the min-max problem. 

For the min-sum problem, if the points in S are given initially all lying on the boundary of C, 
then this case is referred to as the boundary case of the min-sum problem. 

1.2 Applications in Wireless Sensor Networks 

A Wireless Sensor Network (WSN) is composed of a large number of sensors which monitor some 
surrounding environmental phenomenon. Usually, the sensors are densely deployed either inside the 
target phenomenon or are very close to it jl]. Each sensor is equipped with a sensing device with 
limited battery-supplied energy. The sensors process data obtained and forward the data to a base 
station. A typical type of WSN applications is concerned with security and safety systems, such as 
detecting intruders (or movement thereof) around infrastructure facilities and regions. Particularly, 
it is often used to monitor a protected area so as to detect intruders as they penetrate the area or 
as they cross the area border. For example, research efforts have been under way to extend the 
scalability of wireless sensor networks to the monitoring of international borders |10[ [T3] . 

The study of barrier coverage using mobile sensors was originated in [5l [13] and later in [2]. 
Different from the traditional concept of full coverage, it seeks to cover the deployment region by 
guaranteeing that there is no path through the region that can be traversed undetectedly by an 
intruder, i.e., all possible crossing paths through the region are covered by the sensors [2} [5l [T3]. 
Hence, an interesting problem is to reposition the sensors quickly so as to repair the existing security 
hole and thereby detect intruders [2]. Since barrier coverage requires fewer sensors for detecting 
intruders, it gives a good approximation of full area coverage. The planar region on which the 
sensors move is sometimes represented by a circle. Since sensors have limited battery-supplied 
energy, we wish to minimize their movement. Thus, if each sensor is represented as a point, the 
problem is exactly our optimal point movement min-max (the optimization version) or min-sum 
problem. Further, if each sensor has energy A and we want to determine whether this level of 
energy is sufficient to form a barrier coverage, then the problem becomes the decision version of 
the min-max problem. 

1.3 Previous Work and Our Results 

For the min-max problem, Bhattacharya et al. [2] proposed an 0(n^'^) time algorithm for the 
decision version and an 0{n^'^ logn) time algorithm for the optimization version, where the decision 
algorithm is based on some observations and brute force and the optimization algorithm is based 
on parametric search approach [TJ |T5]. Recently, it was claimed in [17] that these two problem 
versions were solvable in 0{n'^'^) time and 0(n^'^ logn) time, respectively. However, it seems that 
the announced algorithms in [17j contain errors (which might be fixed, say, by using the methods 
given in this paper). In this paper, we solve the decision version in 0(relog n) time and the 
optimization version in 0(n log n) time, which significantly improve the previous results. The 
improvements of our algorithms are based on new observations and interesting techniques. 

A by-product of our techniques that is interesting in its own right is an algorithm for dynamically 
maintaining the maximum matchings of circular convex bipartite graphs. Our algorithm handles 
each (online) vertex insertion or deletion on an n- vertex circular convex bipartite graph in 0(log n) 
time. This matches the performance of the best known dynamic matching algorithm for convex 
bipartite graphs [3]. Note that convex bipartite graphs are a subclass of circular convex bipartite 
graphs |14j . To our best knowledge, no dynamic matching algorithm for circular convex bipartite 



graphs was known before. Since dynamically maintaining the maximum matching of a graph is a 
basic problem, our result may find other applications. 

For the min-sum problem, an 0{n'^) time approximation algorithm with approximation ratio 
l+vr was given in [2\. A PTAS approximation algorithm, which has a substantially larger polynomial 
time bound, was also given in [2j. In this paper, we present an 0{'n?) time approximation algorithm 
with approximation ratio 3, which improves the (1 + 7r)-approximation result in [2]. However, 
whether the general min-sum problem is NP-hard is still left open. 

For the boundary case of the min-sum problem, an 0{n^) time (exact) algorithm was given in 
|17j . We show that the time bound of that algorithm can be reduced to 0{n^). 

The rest of this paper is organized as follows. Our algorithm for the decision version of the 
min-max problem is given in Section [2l and our algorithm for the optimization version is presented 
in Section [3l The min-sum problem is discussed in Section [H 

To distinguish from a normal point in the plane, in the following paper we refer to each point 
Ai £ S as a sensor. 

2 The Decision Version of the Min-max Problem 

For simplicity, we assume the radius of the circle C is 1. Denote by dC the boundary of C. Let 
Ac be the maximum distance traveled by the sensors in S in an optimal solution for the min-max 
problem, i.e., Ac = min{maxo<i<n-i{|^i^i|}}. Since the sensors are all in C, Ac < 2. In this 
section, we consider the decision version of the min-max problem on C: Given a value A, determine 
whether Ac < A. We present an 0(nlog n) time algorithm for this problem. 

2.1 An Algorithm Overview 

We first discuss some concepts. A bipartite graph G = {Vi, V2,E) with jVi | = 0{n) and IV2I = 0{n) 
is convex on the vertex set V2 if there is a linear ordering on V2, say, V2 = {wqj f^i) • • • ; ^n-i}, such 
that if any two edges {v, Vj) € E and {v, v^) G E with Vj,Vk € V2, v £ Vi, and j < k, then (v, vi) € E 
for all j < I < k. In other words, for any vertex f G Vi, the subset of vertices in V2 connected to v 
forms an interval on the linear ordering of V2 ■ For any v £Vi, suppose the subset of vertices in V2 
connected to v is {vj,Vj+i, . . . , Vk}; then we denote begin(v, G) = j and end(v, G) = k. Although 
E may have 0{n'^) edges, it can be represented implicitly by specifying begin{v,G) and end{v,G) 
for each v £ Vi. A vertex insertion on G is to insert a vertex v into Vi with an edge interval 
[begin{v, G), end{v, G)] and implicitly connect v to every Vi € V2 with begin{v, G) <i < end{v, G). 
Similarly, a vertex deletion on G is to delete a vertex v from Vi as well as all its adjacent edges. 

A bipartite graph G = {Vi,V2,E) is circular convex on the vertex set V2 if there is a circular 
ordering on V2 such that for each vertex v £ Vi, the subset of vertices in V2 connected to v forms 
a circular-arc interval on that ordering. Precisely, suppose such a clockwise circular ordering of V2 
is Vo,vi, . . . ,Vn-i- For any two edges {v,Vj) G E and {v,Vk) € E with Vj,Vk € V2, v £ Vi, and 
j < k, either (v, vi) G E for all j < I < k, 01 {v, vi) € E for all fc < / < n — 1 and {v, vi) G E for all 
< I < j. For each v £ Vi, suppose the vertices of V2 connected to v are from Vj to Vk clockwise 
on the ordering, then begin{v, G) and end{v, G) are defined to be j and k, respectively. Vertex 
insertions and deletions on G are defined similarly. 

A maximum matching in a convex bipartite graph can be found in 0{n) time [9l W\\ I16j . The 
same time bound holds for a circular convex bipartite graph [2]. Brotal et al. [3] designed a data 
structure for dynamically maintaining the maximum matchings of a convex bipartite graph that can 



support each vertex insertion or deletion in 0(Iog n) amortized time. For circular convex bipartite 
graphs, however, to our best knowledge, we are not aware of any previous work on dynamically 
maintaining their maximum matchings. 

The main idea of our algorithm for the decision version of the min-max problem is as follows. 
First, we model the problem as finding the maximum matchings in a sequence of 0{n) circular con- 
vex bipartite graphs, which is further modeled as dynamically maintaining the maximum matching 
of a circular convex bipartite graph under a sequence of 0{n) vertex insertion and deletion oper- 
ations. Second, we develop an approach for solving the latter problem. Specifically, we show that 
the maximum matching of a circular convex bipartite graph of 0{n) vertices can be dynamically 
maintained in 0(log n) time (in the worst case) for each vertex insertion or deletion. Note that 
this result is of independent interest. 

In the following, we first present the problem modeling and then give our algorithm for dynam- 
ically maintaining the maximum matching of a circular convex bipartite graph. 

2.2 The Problem Modeling 

Recall that in the decision version of the min-max problem, our goal is to determine whether 
Ac < A. Let P be an arbitrary regular n-gon with its vertices Pq,Pi, . . . ,-Pn-i ordered clockwise 
on dC. We first consider the following sub-problem: Determine whether we can move all sensors 
to the vertices of P such that the maximum distance traveled by the sensors is at most A. Let Gp 
be the bipartite graph between the sensors Aq, . . . , j4„_i and the vertices of P, such that a sensor 
Ai is connected to a vertex Pj in Gp if and only if |j4jPj| < A. The next lemma is immediate. 

Lemma 1 The bipartite graph Gp is circular convex. 

Proof: This simply follows from the fact that the boundary of any circle of radius A can intersect 
dC at most twice. □ 

To solve the above sub-problem, it suffices to compute a maximum matching M in the circular 
convex bipartite graph Gp (by using the algorithm in [2]). If M is a perfect matching, then the 
answer to the sub-problem is "yes"; otherwise, the answer is "no". Thus, the sub-problem can be 
solved in 0{n) time (note that the graph Gp can be constructed implicitly in 0{n) time, after 
0(n log n) time preprocessing). If the answer to the sub-problem is "yes", then we say that P is 
feasible with respect to the value A. 

If P is feasible, then clearly Ac < A. If P is not feasible, however, Ac > A does not necessarily 
hold, because P may not be positioned "right" (i.e., P may not be the regular n-gon in an optimal 
solution of the optimization version of the min-max problem). To further decide whether Ac < A, 
our strategy is to rotate P clockwise on dC by an arc distance at most Itx jn. Since the perimeter of 
C is 27r, the arc distance between any two neighboring vertices of P is 27r/n. A simple yet critical 
observation is that Ac < A if and only if during the rotation of P, there is a moment (called a 
feasible moment) at which P becomes feasible with respect to A. Thus, our task is to determine 
whether a feasible moment exists during the rotation of P. 

Consider the graph Gp. For each sensor Ai, denote by E{Ai) = {Pj,Pj^i, . . . , P^} the subset of 
vertices of P connected to Ai in Gp , where the indices of the vertices of P are taken as module by 
n. We assume that E(Ai) does not contain all vertices of P (otherwise, it is trivial). Since the arc 
distance from Pj^i to Pj is 27r/n, during the (clockwise) rotation of P, there must be a moment 
after which Pj-i becomes connected to Ai, and we say that Pj-i is added to E(Ai); similarly, there 
must be a moment after which P^ becomes disconnected to Ai, and we say that Pk is removed 
from E{Ai). Note that these are the moments when the edges of Ai (and thus the graph Gp) are 



changed due to the rotation of P. Also, note that during the rotation, all vertices in E{Ai) \ {Pk} 
remain connected to Ai and all vertices in P \ {E{Ai) U {Pj-i}} remain disconnected to Ai. Hence 
throughout this rotation, there are totally n additions and n removals on the graph Gp. If we sort 
all these additions and removals based on the time moments when they occur, then we obtain a 
sequence of 2n circular convex bipartite graphs, and determining whether there exists a feasible 
moment is equivalent to determining whether there is a graph in this sequence that has a perfect 
matching. With the 0{n) time maximum matching algorithm for circular convex bipartite graphs 
of n vertices in [T3] , a straightforward solution for determining whether there is a feasible moment 
would take O(n^) time. 

To obtain a faster algorithm, we further model the problem as follows. Consider the addition of 
Pj-i to E{Ai). This can be done by first deleting the vertex of Gp corresponding to Ai and then 
inserting a new vertex corresponding to Ai with its edges connecting to the vertices in {Pj_i} U 
E{Ai). The removal of P^ from E{Ai) can be handled similarly. Thus, each addition or removal 
on E(Ai) can be transformed to one vertex deletion and one vertex insertion on Gp. If we sort 
all vertex updates (i.e., insertions and deletions) by the time moments when they occur, then the 
problem of determining whether there is a feasible moment is transformed to determining whether 
there exists a perfect matching in a sequence of vertex updates on the graph Gp. In other words, 
we need to dynamically maintain the maximum matching in a circular convex bipartite graph to 
support a sequence of 2n vertex insertions and 2n vertex deletions. This problem is handled in the 
next subsection. 

2.3 Dynamic Maximum Matching in a Circular Convex Bipartite Graph 

In this subsection, we consider the problem of dynamically maintaining the maximum matching in 
a circular convex bipartite graph to support vertex insertions and deletions. We treat all vertex 
updates in an online fashion. 

Let G = {Vi,V2,E) with l^il = 0{n) and IV2I = 0{n) be a circular convex bipartite graph 
on the vertex set V2, i.e., the vertices of V2 connected to each vertex in Vi form a circular-arc 
interval on the sequence of the vertex indices of V2. Suppose V2 = {uoi^i) • • • ,Vn-i} is ordered 
clockwise. Recall that a vertex insertion on G is to insert a vertex v into Vi with an edge interval 
[begin{v , G) , end{v , G)] such that v is (implicitly) connected to all vertices of V2 from begin{v,G) 
clockwise to end{v, G). A vertex deletion is to delete a vertex v from Vi and all its adjacent edges 
(implicitly). Our task is to design an algorithm for maintaining the maximum matching of G to 
support such update operations (i.e., vertex insertions and deletions) efficiently. Below, we present 
an algorithm with an 0(log n) time per update operation. 

Our approach can be viewed as a combination of the data structure in [3j for dynamically 
maintaining the maximum matching in a convex bipartite graph and the linear time algorithm in 
|14j for computing a maximum matching in a circular convex bipartite graph. We refer to them as 
the BGHK data structure [3] and the LB algorithm [T3], respectively. We first briefly describe the 
BGHK data structure and the LB algorithm. 

The BGHK data structure [3j is a binary tree T, and each node of T maintains a balanced 
binary tree. This data structure can be constructed in 0(?ilog n) time and can support each 
vertex insertion or deletion in 0(log n) amortized time. Consider a vertex insertion, i.e., inserting 
a vertex v into Vi. Let M' (resp., M) be the maximum matching in the graph before (resp., after) 
the insertion. Let \M\ denote the number of matched pairs in M. After the data structure is 
updated (in 0(log n) amortized time), the value \M\ can be reported in 0(1) time and M can 
be reported in 0(|M|) time. We can also determine in 0(1) time whether v is matched in M. 



Further, if another vertex v' G Vi was matched in M' but is not matched in M, then it is easy to 
see that v must be matched in M. When this case occurs, we say that v replaces v' and v' is called 
the replacement, and the data structure is able to report the replacement in 0(1) time. Note that 
as shown in [3], although an update on the graph can cause dramatic changes on the maximum 
matching, the sets of the matched vertices in Vi (and V2) can change by at most one vertex. Thus, 
there is at most one such replacement v'. Similarly, consider deleting a vertex v from Vi. After the 
data structure is updated, the value \M\ can be reported in 0(1) time and M can be reported in 
0(|M|) time. We can also find out whether v was matched in M' in 0(1) time. If a vertex v' G Vi 
was not matched in M' but is matched in M, then it is easy to see that v must be matched in M' . 
When this case occurs, we say v' is the supplement, which can be determined in 0(1) time. 

The LB algorithm ^14j finds a maximum matching in a circular convex bipartite graph G = 
(yi,V2,E) by reducing the problem to two sub-problems of computing the maximum matchings 
in two convex bipartite graphs Gi and G2. Some details are summarized below. For any vertex 

V £ Vi, if begin{v,G) < end{v,G), then v is called a non-boundary vertex. Otherwise, f is a 
boundary vertex; the edges connecting v to Vj,egin{v,G)^''^begin{v,G)+i^ ■ ■ ■ )^n-i in ^2 suce called lower 
edges, and the other edges connecting v are upper edges. Based on the graph G, a convex bipartite 
graph Gi = {Vi,V2,Ei) is defined as follows. Both its vertex sets are the same as those in G. 
For each vertex f G Vi in G, begin{v,Gi) = begin{v,G); if w is a non-boundary vertex, then 
end{v,Gi) = end{v,G), and otherwise end{v,Gi) = n — 1 + end{v,G) (note that this value of 
end{v,Gi) is used only for comparison in the algorithm although there are not so many vertices 
in V2). The LB algorithm has two main steps. The first step is to compute a maximum matching 
in Gi, which can be done in 0{n) time [9l W\\ I16j . Let M{Gi) be the maximum matching of 
Gi. Next, another convex bipartite graph G2 = {Vi^V2-:E2) is defined based on M{Gi) and 
G, as follows. Both its vertex sets are the same as those in G. For each non-boundary vertex 
t; G Vi in G, begin{v,G2) = begin{v,G) and end{v,G2) = end{y,G). For each boundary vertex 

V G Fi in G, there are two cases: If v is matched in M{Gi), then begin{v,G2) = begin{y,G) and 
end{v,G2) = n — 1; otherwise, begin{v,G2) = and end{v,G2) = end(v,G). The second step of 
the LB algorithm is to compute a maximum matching in G2 (in 0{n) time), denoted by M{G2)- 
It was shown in ^14] that M{G2) is also a maximum matching of the original graph G. 

We now discuss our algorithm for dynamically maintaining a maximum matching in the circular 
convex bipartite graph G. As preprocessing, we first run the LB algorithm on G, after which both 
the convex bipartite graphs Gi and G2 of G are available. We then build two BGHK data struc- 
tures for Gi and G2, denoted by T{Gi) and T{G2), respectively, for maintaining their maximum 
matchings. This completes the preprocessing, which takes 0(n log n) time. In the following, we 
discuss how to perform vertex insertions and deletions. 

Consider a vertex insertion, i.e., inserting a vertex v into Vi with the edge interval \begin{v, G), 
end{v, G)]. To perform this insertion, intuitively, we need to update the two BGHK data structures 
T{Gi) and T{G2) in a way that mimics some behavior of the LB algorithm. Specifically, we first 
insert v into the graph Gi by updating T{Gi). Based on the results on Gi (e.g., whether there is a 
replacement) and the behavior of the LB algorithm, we modify G2 by updating T{G2) accordingly. 
In this way, the maximum matching maintained by T{G2) is the maximum matching of G after the 
insertion. The details are given below. 

Let G'l and G'2 be the two graphs that would be produced by running the LB algorithm on G 
with the new vertex v (and its adjacent edges). Let M{Gi),M{G2),M{G'i), and M{G'2) be the 
maximum matchings of Gi,G2,G'i, and G^, respectively. Depending on whether f is a boundary 
vertex, there are two main cases. 



If f is a non-boundary vertex (i.e., begin{v,G) < end{v,G)), then G'l can be obtained by 
inserting v into Gi. Hence we insert v into T{Gi). Depending on whether there is a replace- 
ment, there are two cases. 

— If no replacement, then G2 can be obtained by inserting v into G2. Thus, we simply 
insert v into T{G2) and we are done. 

— Otherwise, let v' be the replacement. So v' was matched in M{Gi) but is not matched 
in M{G'i). Depending on whether v' is a boundary vertex, there are two subcases. 

* If v' is a non-boundary vertex, then again, G'2 can be obtained by inserting v into 
G2. We thus insert v into T{G2) and we are done. 

* If v' is a boundary vertex, then since v' was matched in M(Gi), according to the LB 
algorithm, v' with the edge interval [begin{v' , G),n — 1] is in G2. After the insertion 
of V into Gi, v' is not matched in M{G'i). Thus, according to the LB algorithm, 
G2 can be obtained by deleting v' (with the edge interval [begin{v' , G),n — 1]) from 
G2, inserting v' with the edge interval [0, end(v' , G)] into G2, and finally inserting v 
into G2. 

In summary, for this subcase, we delete v' (with the edge interval [begin{v' , G), n—1]) 
from T{G2) and insert v' with the edge interval [0, end(f', G)] into T(G2). Finally, 
we insert v into T{G2), and we are done. 

If t; is a boundary vertex (i.e., begin{v, G) > end{v, G)), then according to the LB algorithm, 
G'l can be obtained by inserting v with the edge interval [begin{v,G),n — 1 -|- end{v,G)] 
into Gi. Thus we insert v with the edge interval [begin{v, G),n — 1 + end{v, G)] into T{Gi). 
Depending on whether there is a replacement, there are two cases. 

— If no replacement, then depending on whether v is matched in M(G'^), there are two 
subcases. 

* If f is matched, then according to the LB algorithm, G'2 can be obtained by inserting 
V with the edge interval [begin{v, G),n — 1] into G2. Thus, we insert v with the edge 
interval [begin(v , G) , n — 1] into T{G2), and we are done. 

* If w is not matched, then according to the LB algorithm, G'2 can be obtained by 
inserting v with the edge interval [0,end{v,G)] into G'2. We thus insert v with the 
edge interval [0,end{v,G)] into T{G2), and we are done. 

— Otherwise, there is a replacement v'. So v' was matched in M{Gi) but is not matched 
in M{G'i), and v is matched in M{G'i). Depending on whether v' is a boundary vertex, 
there are two subcases. 

* If v' is a non-boundary vertex, then since v is matched in M(G'i), G'2 can be obtained 
by inserting v with the edge interval [begin{v , G) , n — 1] into G2. We thus insert v 
with the edge interval [begin{v , G) , n — 1] into T{G2)- 

* If v' is a boundary vertex, then according to the LB algorithm, G'2 is the graph 
obtained by deleting v' (with the edge interval [begin{v' , G),n—1]) from G2, inserting 
v' with the edge interval [0, end{v' , G)] into G2, and finally inserting v with the edge 
interval [begin{v , G) , n — 1] into G2. 

Thus, we delete v' (with the edge interval [begin{v' , G) , n — 1]) from T(G2), and 
insert v' with the edge interval [0,end{v' ,G)] into T(G2). Finally, we insert v with 
the edge interval [begin{v , G) , n — 1] into T(G2). 

7 



This completes the description of our procedure for handling a vertex insertion. 

Next, consider a vertex deletion, i.e., deleting a vertex v from Vi of G. Our procedure for this 
operation proceeds in a manner symmetric to the insertion procedure, and we briefly discuss it 
below. Define the two graphs G'l and G2 similarly as above. 

• If f is a non-boundary vertex, then we delete v from T(Gi). If no supplement, then we delete 
V from T{G2) and we are done. Otherwise, let v' be the supplement. So v' was not matched 
in M[Gi) but is matched in M{G'i). Depending on whether v' is a boundary vertex, there 
are two cases. 

— If v' is a non-boundary vertex, then we delete v from T{G2) and we are done. 

— If v' is a boundary vertex, then we delete v' with the edge interval [0,end(v' ,G)] from 
T{G2) and insert v' with the edge interval [begin{v' , G),n — 1] into T{G2)- Finally, delete 

V from T(G2), and we are done. 

• If w is a boundary vertex, then we delete v (with the edge interval [begin{v , G) , n — 1 + 
end{v,G)]) from T{Gi). Depending on whether there is a supplement, there are two cases. 

— If no supplement, then depending on whether v was matched in M(Gi), there are two 
subcases. If v was matched, then we delete v (with the edge interval [begin{v, G),n — 1]) 
from T(G2); otherwise, we delete v (with the edge interval [0,end{v,G)]) from T{G2)- 

— Otherwise, let v' be the supplement. So v' was not matched in M(Gi) but is matched 
in M{Gi), and v was matched in M{Gi). Since v was matched in M{Gi), according to 
the LB algorithm, G2 contains v with the edge interval [begin{v , G) , n — 1]. If v' is a 
non-boundary vertex, then we delete v (with the edge interval [begin{v, G),n — 1]) from 
T{G2) and we are done. Otherwise, since v' was not matched in M(Gi), according to the 
LB algorithm, G2 contains v' with the edge interval [0,end{v' ,G)]; since v' is matched 
in M{G'i), according to the LB algorithm, G2 should contain v' with the edge interval 
[begin{v' , G),n — 1]. Therefore, we delete v' (with the edge interval [0, end{v' , G)]) from 
T{G2), insert v' with the edge interval [begin{v' , G),n — 1] into T{G2), and finally delete 

V (with the edge interval [begin(v , G) , n — 1]) from T{G2)- 

This completes the description of our vertex deletion procedure. 

As shown in Subsection 12.21 the decision version of the min-max problem can be transformed 
to the problem of dynamically maintaining the maximum matching in a circular convex bipar- 
tite graph subject to a sequence of vertex insertions and deletions. Hence, the correctness of our 
algorithm for the decision version hinges on the correctness of our dynamic maximum matching 
algorithm for circular convex bipartite graphs. Yet, the correctness of our (online) dynamic max- 
imum matching algorithm for circular convex bipartite graphs can be seen quite easily. This is 
because our procedures for performing vertex insertions and deletions are both based on the fact 
that they simply mimic the behavior of the LB algorithm (while implementing their processing by 
the means of the BGHK data structures). 

For the running time of our algorithm, each update operation involves at most two vertex 
insertions and two vertex deletions on T{G\) and T{G2), each of which takes 0(log n) amortized 
time [3]; thus, it takes 0(log n) amortized time in total. Actually, the BGHK data structure in 
[3j supports vertex insertions and deletions not only on Vi but also on ¥2- Inserting vertices on V2 
may make the tree unbalanced, and that is why its running time is amortized. However, if vertices 





(a) (b) 

Figure 1: (a) The points Xi and Yi on dC for Ai; (b) \AiA'-\ = \AjA'-\. 



are inserted only on Vi, then the tree will never become unbalanced and thus each update takes 
0(log n) time in the worst case. In our problem formulation, the vertex updates indeed are only 
on Vi. Denote by M{G) the maximum matching in G. We then have the following result. 

Theorem 1 A data structure on a circular convex bipartite graph G = (Vi,V2,-E') can he built in 
0(n log n) time for maintaining its maximum matching M{G) so that each online vertex insertion 
or deletion on Vi can be done in 0(log n) time in the worst case. After each update operation, 
\M{G)\ can be reported in 0(1) time and M{G) can be reported in 0{\M{G)\) time. 

Since the decision version of the min-max problem has been reduced to dynamically maintaining 
the maximum matching in a circular convex bipartite graph under a sequence of 2n vertex insertions 
and 2n vertex decisions, we solve the dynamic maximum problem as follows. After each update 
operation, we check whether |M((jr)| = n, and if this is true, then we report Ac < A and halt the 
algorithm. If all An updates have been processed but it is always |M(G)| < n, then we report 
Ac > A. Based on Theorem [H we have the result below. 

Theorem 2 Given a value \, we can determine whether Xc < X in 0(nlog n) time for the decision 
version of the min-max problem. 

3 The Optimization Version of the Min-max Problem 

In this section, we consider the optimization version of the min-max problem, and present an 
0(n log n) time algorithm for it. The main task is to compute the value Ac. 

Let o be the center of G. For simplicity of discussion, we assume that no sensor lies at o. Denote 
by Xi and Yi the two points on dG which are closest and farthest to each sensor Ai, respectively. 
Clearly, Xi and Yi are the two intersection points of dC with the line passing through Ai and the 
center o of C (see Figure 1(a)). The lemma below has been proved in [IT], and for self-containment 
of this paper, we include that proof in Appendix lAl 

Lemma 2 JJ7| / Suppose an optimal solution for the min-max optimization problem is achieved with 
Ac = l^j^il for some i G {0, . . . , n — 1}. Then either A'- is the point Xi, or there is another sensor 



^j (J 7^ i) such that Xc 



\A- A' \ 



also holds. In the latter case, any slight rotation of the regular 



n-gon that achieves Xc in either direction causes the value of Xc to increase (i.e., it makes one of 
the two distances \AiA'-\ and \AjA'A increase and the other one decrease). 

The points on dC satisfying the conditions specified in Lemma [2] may be considered as those 
defining candidate values for Ac, i.e., they can be considered as some vertices of possible regular 
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n-gons on dC in an optimal solution. The points Xh of all sensors Ah (0 < /i < n — 1) can be 
easily determined. Define Di = ^h=o{\^hXh\}, which can be computed in 0{n) time. But, the 
challenging task is to handle all the pairs {Ai,Aj) [i ^ j) such that the distance from Ai to a 
vertex of a regular n-gon is equal to the distance from Aj to another vertex of that n-gon and a 
slight rotation of the n-gon in either direction monotonically increases one of these two distances 
but decreases the other. We refer to such distances as the critical equal distances. Denote by 
D2 the set of all critical equal distances. Let D = Di Li D2. By Lemma [21 Ac € D. Thus, if 
D is somehow available, then Ac can be determined by using our algorithm in Theorem [5] in a 
binary search process. Since Di is readily available, the key is to deal with D2 efficiently. An easy 
observation is maxo</i<n-i |A/iX/i| < Ac- We can use the algorithm in Theorem[2]to check whether 
Ac < maxo</i<„_i |A/iX/t|, after which we know whether Ac = maxo</i<n-i |A/iX/i|. Below, we 
assume maxo</t<„_i |A/iX/i| < Ac (otherwise, we are done). Thus, we only need to focus on finding 
Ac from the set D2 ■ 

It has been shown in [17] that I-D2I = O(n^). Of course, our goal is to avoid an 0{rfi) time 
solution. To do so, first we determine a subset D'2 of D2 such that Ac G -Dg but with jZ^gl = O(n^). 
Furthermore, we do not compute D2 explicitly. Specifically, our idea is as follows. We show that 
the elements of D2 are the y-coordinates of a subset of intersection points among a set F of 0{n) 
functional curves in the plane such that each curve is x-monotone and any two such curves intersect 
in at most one point at which the two curves cross each other. (Such a set of curves is sometimes 
referred to as pseudolines in the literature.) Let Ap be the arrangement of F and \Af\ be the 
number of vertices of Ap- Without computing Ap explicitly, we will generalize the techniques in 
[8] to compute the A:-th highest vertex of Ap for any integer k with 1 < k < \Ap\ in 0(n log n) 
time. Consequently, with Theorem [21 the value Ac can be computed in 0(n log n) time. The 
details are given below. 

Let P be an arbitrary regular n-gon with its vertices Pq, Pi, . . ., Pn~i clockwise on dC. Suppose 
the distances of all the pairs between a sensor and a vertex of P are di < d2 < ■ ■ ■ < (i„2 in sorted 
order. Let do = 0. Clearly, do < Ac < d„2 (the case of Ac = is trivial). Hence, there exists an 
integer k with Q <k <v? such that Ac € {dk,dk+i\- One can find dk and dfc+i by first computing 
all these n^ distances explicitly and then utilizing our algorithm in Theorem [2] in a binary search 
process. But that would take r2(n^) time. In the following lemma, we give a faster procedure 
without having to compute these ii? distances explicitly. 

Lemma 3 The two distances dk and dk+i can he obtained in 0(n log n) tinne. 

Proof: We apply a technique, called binary search in sorted arrays [6], as follows. Given M arrays 
Ai, 1 < i < M, each containing 0{N) elements in sorted order, the task is to find a certain element 
5 ^ A = U^^Aj. Further, assume that there is a "black-box" decision procedure 11 available, such 
that given any value a, 11 reports a < 6 01 a> 6 va 0(T) time. An algorithm is given in [6] to find 
the sought element 6 in A = ufl^Ai in 0{(M + T) \og{NM)) time. We use this technique to find 
dk and dfc+ii as follows. 

Consider a sensor A^. Let S{Ai) be the set of distances between Ai and all vertices of P. In 
O(logn) time, we can implicitly partition S{Ai) into two sorted arrays in the following way. By 
binary search, we can determine an index j such that Xi lies on the arc of dC from Pj to Pj+i 
clockwise (the indices are taken as module by n). Recall that Xi is the point on dC closest to A^. 
If a vertex of P is on Xi , then define j to be the index of that vertex. Similarly, we can determine 
an index h such that Yi (i.e., the farthest point on dC to Ai) lies on the arc from P^ to Ph+i 
clockwise. If a vertex of P is on Yi, then define h to be the index of that vertex. Both j and 
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h can be determined in O(logn) time, after which we imphcitly partition S{Ai) into two sorted 
arrays: One array consists of all distances from Ai to Pj,Pj-i, . . . , Ph+i, and the other consists of 
all distances from Ai to Pj+i, Pj+2, ■ ■ ■ ,Ph (again, all indices are taken as module by n). Note that 
both these arrays are sorted increasingly and each element in them can be obtained in 0(1) time 
by using its index in the corresponding array. 

Thus, we obtain 2n sorted arrays (represented implicitly) for all n sensors in 0(n log n) time, 
and each array has no more than n elements. Therefore, by using the technique of binary search in 
sorted arrays, with our algorithm in Theorem [2] as the black-box decision procedure, both d^. and 
dk+i can be found in 0(n log n) time. The lemma thus follows. □ 

By applying Lemma O we have Xc € {dk,dk+i]- Below, for simplicity of discussion, we as- 
sume Ac ^ dk+i- Thus Ac G {dk,dk+i)- Since maxo<h<n-i\AhXh\ < Xc, we redefine dk := 
max{(ifc,maxo</i<n-i |Ah-^/i|}- We still have Ac € {dk, dk+i)- Let L'2 be the set of all critical equal 
distances in the range (dfc,(ifc+i)- Then Ac € -Dj. We show below that ID2I = 0{n'^) and Ac can 
be found in O(nlog^n) time without computing Dg explicitly. 

Suppose we rotate the regular n-gon P = {Pq, Pi, ... , P„_i) on dC clockwise by an arc distance 
27r/n (this is the arc distance between any two adjacent vertices of P). Let Ai{Pfi{t)) denote the 
distance function from a sensor Ai to a vertex Ph of P with the time parameter t during the 
rotation. Clearly, the function Ai{Ph{t)) increases or decreases monotonically, unless the interval 
of dC in which Ph moves contains the point Xj or Yf, if that interval contains Xj or Yi, then we 
can further divide the interval into two sub-intervals at Xj or 1^, such that Ai{Ph{t)) is monotone 
in each sub- interval. The functions Ai{Ph(t)), for all P^s of P, can thus be put into two sets Sn 
and Si2 such that all functions in Sn monotonically increase and all functions in Si2 monotonically 
decrease. Let m = \Sii\. Then m < n. Denote hy d^ < d^ < ■ ■ ■ < d\j^ the sorted sequence of the 
initial values of the functions in Sn. Also, let dg = and (i^+i = 2 (recall that the radius of C 
is 1). It is easy to see that the range {dk,dk+i) obtained in Lemma[3]is contained in [(i*-,(iV -^] for 
some < j < m. The same discussion can be made for the distance functions in the set Si2 as well. 

Since we rotate P by only an arc distance 27r/n, during the rotation of P, each sensor Ai can 
have at most two distance functions (i.e., one decreasing and one increasing) whose values may 
vary in the range {dk,dk+i). We can easily identify these at most 2n distance functions for the n 
sensors in 0(n log n) time. Denote by F' the set of all such distance functions. Clearly, all critical 
equal distances in the range (d^, d^+i) can be generated by the functions in F' during the rotation 
of P. Because every such distance function either increases or decreases monotonically during the 
rotation of P, each pair of one increasing function and one decreasing function can generate at most 
one critical equal distance during the rotation. (Note that by Lemma El a critical equal distance 
cannot be generated by two increasing functions or two decreasing functions.) Since \F'\ < 2n, the 
total number of critical equal distances in {dk,dk+i) is bounded by 0{n'^), i.e., \D2\ = 0{v?). For 
convenience of discussion, since we are concerned only with the critical equal distances in [dk, dk+i), 
for each function in F' , we restrict it to the range {dk,dk+i) only. 

Let the time t be the x-coordinate and the function values be the y-coordinates of the plane. 
Then each function in F' defines a curve segment that lies in the strip of the plane between the 
two horizontal lines y = d^ and y = dfc+i- We refer to a function in F' and its curve segment 
interchangeably, i.e., F' is also a set of curve segments. Clearly, a critical equal distance generated 
by an increasing function and a decreasing function is the y-coordinate of the intersection point of 
the two corresponding curve segments. Note that every function in F' has a simple mathematical 
description. Below, we simply assume that each function in F' is of 0(1) complexity. Thus, many 
operations on them can each be performed in 0(1) time, e.g., computing the intersection of a 
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decreasing function and an increasing function. 

The set D2 can be computed explicitly in 0{n^) time, after which Ac can be easily found by 
binary search. Below, we develop a faster solution without computing D2 explicitly, by utilizing 
the property that each element of D2 is the y-coordinate of the intersection point of a decreasing 
function and an increasing function in F' and generalizing the techniques in [8]. 

A slope selection algorithm for a set of points in the plane was given in [8]. We will extend this 
approach to solve our problem. The following lemma is needed. 

Lemma 4 For any two increasing (resp., decreasing) functions in F' , if the curve segments defined 
by them are not identical to each other, then the two curve segments intersect in at most one point 
and they cross each other at their intersection point (if any). 

Proof: We only prove the decreasing case. The increasing case can be proved similarly. Let 
Ai[Pa{t)) and Aj{Pb{t)) be two decreasing curves in F\ where Ai{Pa{t)) (resp., Aj{Pb{t))) is the 
distance function between the sensor Ai (resp., Aj) and the vertex Pa (resp., P},) of the regular 
n-gon P, and the two curve segments defined by Ai{Pa{t)) and Aj[Pi,{t)) are not the same. Since 
each sensor has at most one decreasing function in F' , we have Ai 7^ Aj. We assume that during 
the (clockwise) rotation of P, Ai{Pa{t)) = Aj(Pi,(t)) at the moment t = ti and ti is the first such 
moment. Below, we prove that Ai{Pa{t)) = Aj{Pb(t)) cannot happen again for any t > ti in the 
rotation. There are two cases: Pa = Pb and Pa j^ Pb- 

For any two points p and q, let l{p,q) denote the line passing through the two points and pq 
denote the line segment with endpoints p and q whose length is \pq\. Recall that o is the center of 
the circle C. Let Pa{ti) and ^^(ti) be the positions of Pa and Pb at the moment ti, respectively. 

• Pa = Pb- Clearly, Paiti) = Pbih). Let / be the perpendicular bisector of the line segment 
AiAj. At the moment ti, since |^jPa(ti)| = \AjPa{ti)\, Pa{ti) is at one of the two intersection 
points of I and dC. Further, since Ai{Pa{t)) is a decreasing function, Pa{ti) must be on the 
right side of the line l{Ai,o) if we walk from Ai to o (see Fig. EJ^a)). Similarly, Pa{ti) must 
be on the right side of the line l{Aj,o) (going Aj to o). Let z be the other intersection point 
of I and dC. It is easy to see that z is on the left side of either the line l{Ai,o) or the line 
l{Aj,o). Note that d^ > maxo</i<„„i |A/iXft|. Thus, d^ > |^j^i| and dk > |^jXj|. During 
the rotation of P, since both Ai{Pa{t)) and Aj{P{j{t)) are always larger than d^, Pa{t) cannot 
pass any of Xi and Xj, and thus Pa{t) cannot arrive to the position z during the rotation. 
Hence, Ai{Pa{t)) = Aj{Pb{t)) cannot happen again after ti. 

Further, recall that ti is the first moment from the beginning of the rotation with Ai{Pa{t)) = 
Aj{Piy{t)). Without loss of generality, we assume Ai{Pa{t)) < Aj{Ph{t)) for any time t < ti 
(as the example shown in Fig.[2ja)). It is easy to see that Ai{Pa{t)) > Aj{Ph(t)) for any time 
t > ti, which implies that the two functions cross each other at their intersection point. 

• Pa ^ Pfy. At the moment ti, we have |AjPa(ti)| = |AjPfe(ii)|. Assume to the contrary that 
at some moment ^2 > ti, we also have Ai(Pa{t2)) = Aj{Pi,{t2))- Suppose at the moment t2, 
Paih) is at the position P^ and Pb{t2) is at the position P^' (see Fig. [5]|^b)). Then |AjP^| = 
|^jP^'|. Since Pa and Pb are rotated simultaneously, the arc distance from Pa(ti) to P^ is 
equal to the arc distance from Pbiti) to P^', and thus \Pa{ti)Pa\ = \Pb{ti)Pl\. Consider the 
two triangles APb{ti)AjPl^ and APa(ti)AiPa (shown with red solid segments in Fig. E^b)). 
Since \AiPa{ti)\ = \AjPb{ti)\, \AiP^\ = \AjPl\, and |Pa(ii)P^| = \Pb{ti)Pl,\, APb{ti)AjPl is 
congruent to APa{ti)AiP^. Thus, the two angles AAiPa{ti)P'a = ZAjPb{ti)Pl. Further, it is 
easy to see ZoPaiti)Pa = ^oPb{t\)P^. Consequently, we have ZoPa(ti)Aj = AoP},{t\)Aj. 
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Figure 2: Illustrating the proof of Lemma [H (a) Pa = Pf, (b) Pa -^ Pb- 

But, if ZoPa{ti)Ai = ZoPf,{ti)Aj, then we can show that the two functions Ai{Pa{t)) and 
Aj{Pii{t)) define exactly the same curve segment. The proof is nothing but the inverse 
of the above argument. Specificahy, consider any time moment t^ > ti before the end 
of the rotation. Suppose at the moment t^, Pa is at the position P^ and Pb is at the 
position Pi'. Since ZoPa{ti)Ai = ZoPbiti)Aj and ZoPa{ti)P'^ = ZoPbiti)Pl', we have 
ZAiPa{ti)P'^ = ZAjPbiti)Pl'. Further, since \AiPa{ti)\ = \AjPbiti)\ and |Pa(ti)P^'| = 
\Pb{ti)Pl'\, APb{ti)A,Pl'iscongvnenttoAPa{ti)AiP;^. Thus, \AiP^\ = \AjPl'\,le., Ai{Pa{t)) 
Aj{Pb{t)) at any time t = t^ > ti. Similarly, we can also show that at any time moment 
ts < ti, Ai{Pa{t3)) = Aj{Pb{t3)). Hence, Ai{Pa{t)) and Aj{Pb{t)) define exactly the same 
curve segment. But this contradicts with the fact that the curve segments defined by these 
two functions are not the same. This implies that Ai{Pa{t)) = Aj{Pb{t)) cannot happen again 
at any moment t > ti. 

Further, without loss of generality, we assume Ai(Pa(t)) < Aj{Pb{t)) for any time t < ti (as 
the example shown in Fig. E^b)). We then show that Ai{Pa{ts)) > Aj{Pb{t3)) for any time 
is > ti, which means that the two functions cross each other at their intersection point. We 
briefly discuss this. Again, suppose at the moment is. Pa is at the position P^' and Pb is at the 
position Pi'. First, since Ai{Pa{t)) < Aj{Pb{t)) for any time t < ti, it must be \oAj\ > \oAi\ 
(this can be proved by similar techniques as above and we omit the details). Consider the 
two triangles AoAiPa{ti) and AoAjPb{ti) (at the moment ti). Since |AjPa(^i)| = l^jA(ii)l) 
|o-Pa(ii)| = lo-PfoC^i)!, and \oAj\ > \oAi\, we have ZoPb{ti)Aj > ZoPa{ti)Ai, which further 
implies ZAjPbiti)Pl' < AiPa{ti)P'^. Consider the triangles APb{ti)AjPll and APa{ti)AiP'^. 
Due to \AiPa{ti)\ = \AjPb{tx)\, \Pa{ti)P'^\ = \Pbiti)Pl'\, and ZA,Pbiti)Pl' < AiPa{ti)P'^, it 
must be |^jP^'| < \AiP'J^\. In other words, Ai{Pa{t)) > Aj{Pb{t)) at any time t = ts > ti. 

The lemma thus follows. □ 

We further extend every curve segment in F' into an x-monotone curve, as follows. For each 
increasing (resp., decreasing) curve segment, we extend it by attaching two half-lines with slope 1 
(resp., —1) at the two endpoints of that curve segment, respectively, such that the resulting new 
curve is still monotonically increasing (resp., decreasing). Denote the resulting new curve set by 
F. Obviously, an increasing curve and a decreasing curve in F intersect once and they cross each 
other at their intersection point. For any two different increasing (resp., decreasing) curves in F, 
by Lemma U] and the way we extend the corresponding curve segments, they can intersect in at 
most one point and cross each other at their intersection point (if any). In other words, F can be 
viewed as a set of pseudolines. Let Af be the arrangement of F. Observe that the elements in D'2 
are the y-coordinates of a subset of the vertices of Af- Since Xc & D'2, Xc is the y-coordinate of a 
vertex of Af- Denote by \Af\ the number of vertices in Af- Of course, we do not want to compute 
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the vertices of Ap explicitly. By generalizing some techniques in [8] , we have the following lemma. 

Lemma 5 The value \Af\ can be computed in 0(n log n) time. Given an integer k with 1 < k < 
\Af\, the k-th highest vertex of Ap can he found in 0(n log n) time. 

Proof: First of all, because every function in F' is of 0(1) complexity, we can determine in 0(1) 
time whether a curve segment in F' intersects a given horizontal line, and if "yes" , then compute the 
intersection. Thus, for every curve in F, we can also compute its intersection with any horizontal 
line in 0(1) time. Let A^ = |F| < 2n. 

Recall that the curve segments in F' are all in the horizontal strip between y = dk and y = 
dk+i- Thus, all vertices of Ap above the horizontal line y = d^+i are intersections of the newly 
attached half-lines. We can easily determine the highest vertex of Ap in O(nlogn) time, e.g., 
by using the approach in p]. Let / be a horizontal line higher than the highest vertex. Denote 
by /i , /2 , • • • , /tv the sequence of the curves of F sorted in increasing order of the x-coordinates 
of their intersections with I. Similarly, we can determine the lowest vertex of Ap in O(nlogn) 
time. Let f-n{i)-: fTT(2)i ■ ■ ■ ^ fniN) be the sequence of the curves of F sorted in increasing order of 
the x-coordinates of their intersections with a horizontal line below the lowest vertex of Ap. Since 
the curves in F can be viewed as a set of pseudolines, as in [S], the number of inversions in the 
permutation vr, which can be computed in O(nlogn) time, is equal to \Ap\. In summary, we can 
compute \Ap\ in 0(n log n) time. 

To compute the A:-th highest vertex of Ap, we choose to generalize the 0{n log n) time algorithm 
in [8]. Let L be a set of n lines in the plane and Al be the arrangement of L. An 0(n log n) 
time algorithm was given in [8j for computing the k-th highest vertex of Al (1 < ^ < I^lI) in 
0(n log n) time based on parametric search [3 [15]. The main property used in the algorithm [8j is 
the following one. Denote hy li,l2, ■ ■ ■ ,ln the sequence of lines in L sorted in increasing order of their 
intersections with a horizontal line above the highest vertex of Al- Given any horizontal line I', 
let Zjr(i) , Z,r{2) ) • • • ) C(n) be the sequence of lines of L sorted in increasing order of their intersections 
with /'. Then, the number of vertices of Al above I' is equal to the number of inversions in the 
permutation vr. 

In our problem, since any two curves in F can intersect each other in at most one point and 
they cross each other at their intersection point, the above property still holds for Ap. Thus, the 
0(n log n) time algorithm in [8] is applicable to our problem. Therefore, we can find the A;-th 
highest vertex of Ap in 0(n log n) time, and the lemma follows. □ 

A remark: An optimal 0(n log n) time algorithm was also given in [8] (and in |12] ) for finding 
the A;-th highest vertex of Al. However, these algorithms are overly complicated. Although we 
think that the 0(n log n) time approach in [8] may be made work for our problem, it does not 
benefit our overall solution for the optimization version of the min-sum problem because its total 
time is dominated by other parts of the algorithm. Hence, the much simpler 0(n log n) time 
solution (for finding the k-th highest vertex of Ap) suffices for our purpose. 

Recall that Ac is the y-coordinate of a vertex of Ap. Our algorithm for computing Ac then 
works as follows. First, compute |^_f|. Next, find the (|^i;'|/2)-th highest vertex of Ap, and denote 
its y-coordinate by A^. Determine whether Ac < Am by the algorithm in Theorem [21 after which 
one half of the vertices of Ap can be pruned away. We apply the above procedure recursively on 
the remaining vertices of Ap, until Ac is found. Since there are O(logn) recursive calls to this 
procedure, each of which takes 0(nlog n), the total time for computing Ac is 0(nlog n). 

Theorem 3 The min-max optimization problem is solvable in 0(n log n) time. 
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4 The Min-sum Problem 

In this section, we present our new algorithms for the min-sum problem. We show that the boundary 
case of this problem is solvable in 0{n'^) time, which improves the 0{n'^) time result in [T7|. We 
also give an 0{in?) time approximation algorithm with approximation ratio 3, which improves the 
(1 + 7r)-approximation 0{'n?) time algorithm in ^. 

For the boundary case, the 0{n^) time algorithm in [T7] uses the 0{n^) time Hungarian algo- 
rithm to compute a minimum weight perfect matching in a complete bipartite graph. However, the 
graph for this case is very special in the sense that all its vertices lie on the boundary of a circle. 
By using the result in [1], we can actually find a minimum weight perfect matching in such a graph 
in 0{n) time. Therefore, if we follow the algorithmic scheme in [T7j but replace the Hungarian al- 
gorithm by the algorithm in [Ij, the boundary case can be solved in 0{'n?') time. For completeness, 
more details on the proof of the following theorem are given in Appendix [Bl 

Theorem 4 The boundary case of the min-sum problem can be solved in 0{n'^) time. 

Next, we discuss our approximation algorithm for the general min-sum problem. 

Let ^0) • • • 5^n-i be the sensors in C. Our approximation algorithm works as follows. (1) For 
each sensor Ai, i = 0, . . . ,n — 1, compute the point Xi on dC that is closest to Ai. (2) By using 
the algorithm in Theorem HI solve the following min-sum boundary case problem: Viewing the n 
points Xq,Xi, . . . ,Xn-i as pseudo- sensors (which all lie on dC), find n points on dC as the goal 
positions for the pseudo-sensors such that the sum of the distances traveled by all n pseudo-sensors 
is minimized. Let X[ be the goal position for each Xj (0 < i < n — 1) in the optimal solution thus 
obtained. We then let X[ be the goal position for each sensor Ai, < i < n — 1, for our original 
min-sum problem. This completes the description of our approximation algorithm. 

Clearly, with Theorem HI the time complexity of the above approximation algorithm is O(n^). 
The lemma below shows that the approximation ratio of this algorithm is 3. 

Lemma 6 The approximation ratio of our approximation algorithm is 3. 

Proof: Let A = Yl'^Za l^i^il- Let Al,A\, . . . ,^*_i be the goal positions of ah sensors (i.e., A* is 
the goal position for each sensor 74j, < i < n— 1) in an optimal solution for the min-sum problem. 
Let Ac = YhZq \A.iA*\. Our task is to prove A < 3 • Ac- 

First, YJlZo \XiX'^ < YJIZq \XiA*\, and \AiX,\ < \A,A*\ holds for each < i < n - 1. Then, 

n— 1 n— 1 



A = y^ \AiXl\ < y^ilAiXil + \XiX-\) {triangle inequality) 

j=0 i=0 

n— 1 n— 1 n— 1 n— 1 

= Y. 1^*^*1 + Yl i^^^^'i ^ Yl i^^^^i + Y. 1^*^**1 

j=0 1=0 i=0 j=0 

n— 1 n— 1 

< 2- > |^jXj| -|- > l^i^i'l {triangle inequality) 

i=Q i=0 

n-l 

< 3-^|Ai4*| = 3-Ac. 



i=0 

The lemma thus follows. □ 

Hence, we conclude with the following result. 
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Theorem 5 There exists an 0{n?) time approximation algorithm for the min-sum problem with 
approximation ratio 3. 
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Appendix 

A The Proof of Lemma [2] 

Lemma [2] [T7] Suppose an optimal solution for the min-max optimization problem is achieved with 
Xc = \AiA'-\ for some i G {0, . . . , n — 1}. Then either A'- is the point Xi, or there is another sensor 
^j (J 7^ ^j such that Xc = \AjA'A also holds. In the latter case, any slight rotation of the regular 
n-gon that achieves Xc in either direction causes the value of Xc to increase (i.e., it makes one of 
the two distances \AiA'j\ and \AjA'-\ increase and the other one decrease). 

Proof: First assume that in an optimal solution, the sensor Ai is the only one satisfying Ac = 
l^jA^I, but A[ is not the point Xi. Thus |AiXj| < |Aj^^| holds (see Figure 1(a)). Then, we 
rotate the regular n-gon that achieves Ac by moving the vertex A[ towards Xj, with a very small 
distance e. Clearly, the distance function between Ai and A[ decreases monotonically during this 
rotation of that n-gon. Denote by Aq, A'l^ . . ., A'l^_^ the new positions of the sensors after the 
rotation stops. Since e is arbitrarily small and Ai is the only sensor satisfying Ac = |AjA^|, we 
have l^jA"! > |Afc^^| for all k ^ i\ moreover, |^iA"| < |AjA^| holds. But, this contradicts with 
the assumption that Ac = l^i^il gives an optimal solution to the min-max optimization problem. 
Suppose now there exists another sensor Aj such that the optimal value Ac = l^j^jl (j 7^ i) also 
holds (see Figure 1(b)). A slight rotation of the regular n-gon that achieves Ac in either direction 
cannot make WjA'-'l < l^i^il and lA,^''! < lAo^'l both occur, where A'J and A'- are the new 
positions of A'^ and A'; after the rotation stops (otherwise, it would contradict with the assumption 
that Ac = l^i^il = l^j^yl gives an optimal solution). Hence, the rotation of the regular n-gon 
that achieves Ac increases one of the two distances l^i^^l and l^djA'!, while decreasing the other. 
The proof is thus complete. □ 

B The Proof of Theorem [4] 

Recall that in the boundary case of the min-sum problem, all sensors are on the boundary dC 
of C. Let Aq, Ai, ..., A^-i denote the initial positions of the n sensors on dC, and A'q, A'^, 
. . ., A'^_i denote their goal positions on dC that form a regular n-gon. Denote by Ac the sum 
of the distances traveled by all n sensors in an optimal solution of the min-sum problem, i.e., 
Ac = ininX^i^o l^j^il- '^^^ following lemma has been proved in [17]. 

Lemma 7 JJ7| / There exists an optimal solution for the boundary case of the min-sum problem 
with the following property: There exists a sensor Ai which does not move, i.e., Ai = A[. 

Based on Lemma[7l the boundary case can be solved as follows. For each sensor Ai, < i < n— 1, 
let P{Ai) be the regular n-gon on dC such that Ai is one of its vertices. Denote by Hi the complete 
bipartite graph between the set of all sensors and the set of all vertices of P{Ai) such that the 
weight of an edge connecting a sensor and a vertex of P{Ai) is defined as their Euclidean distance. 
We compute a minimum weight perfect matching Mi in Hi, for each < i < n — 1, and the one that 
gives the minimum weight defines an optimal solution for our original problem. Here, the weight 
of a perfect matching is the sum of all edge weights of the matching. 

The running time of the above algorithm is dominated by the step of computing the minimum 
weight perfect matchings in the graphs Hi. The algorithm in [TT] uses the 0{n^) time Hungarian 
algorithm for computing such matchings in the graphs Hi. 
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Let H he a complete bipartite graph with two vertex sets of cardinaHties ni and n2, respectively, 
such that all its vertices lie on the boundary of a circle and each edge weight is the Euclidean distance 
between two such vertices (the edges are represented implicitly). A maximum cardinality matching 
of H consists of min{ni, 712} edges. An algorithm was given in [4J for computing a minimum weight 
maximum cardinality matching in H in 0(ni + n2) time (i.e., the total sum of edge weights in the 
output maximum cardinality matching is as small as possible). 

Since in our algorithm, all vertices of every complete bipartite graph Hi lie on dC, the linear 
time algorithm in [2] can be applied to compute a minimum weight maximum cardinality matching 
of Hi in 0{n) time, for < i < n — 1. (Note that a maximum cardinality matching in the graph 
Hi is a perfect matching, and vice versa.) Consequently, the total running time of our algorithm is 
O(n^). Theorem H] thus follows. 
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